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PACS 97 . 10 . Kc - Stellar rotation 

PACS 97.10.Yp - Star counts, distribution, and statistics 

PACS 05 . 90 . +m - Other topics in statistical physics, thermodynamics, and nonlinear dynamical 
systems 

Abstract. - In the present study, we analyze the radial velocity distribution as a function of 
different stellar parameters such as stellar age, mass, rotational velocity and distance to the Sun 
for a sample of 6781 single low-mass field dwarf stars, located in the solar neighborhood. We 
show that the radial velocity distributions are best fitted by q-Gaussians that arise within the 
Tsallis nonextensive statistics. The obtained distributions cannot be described by the standard 
Gaussian that emerges within Boltzmann-Gibbs (B-G) statistical mechanics. The results point 
to the existence of a hierarchical structure in phase space, in contrast to the uniformly occupied 
phase space of B-G statistical mechanics, driven by the g-Central Limit Theorem, consistent with 
nonextensive statistical mechanics. 



Introduction. — The behavior of the distribution of 
different stellar physical parameters, in particular rota- 
tional and radial velocities, seems to be better explained 
on the basis of a Tsallis maximum entropy distribution 
function [I] than by standard Gaussian and Maxwell- 
Boltzmann distributions or analytical functions (Soares et 
al. [5] and Carvalho et al. J3--5J. Indeed, the nonextensive 
generalization of Boltzmann-Gibbs statistics, presented by 
Tsallis PQ, appears now as a powerful parametrization of 
the statistical mechanics of out-of-equilibrium states in 
open systems as is the case of stars, planetary and stellar 
systems. 

The Tsallis statistics, is based on q— exponential and q- 
logarithm functions defined by, 
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from where, associated with the g-statistics emerges the 
entropy S q [B] 
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Fig. 1: The spatial distribuition of the present working sample 
(3) of dwarf stars in the solar neighbourhood. The Sun is centered 
in point (0,0,0). 
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where the g-Gaussian recovers the usual one at q = 
1. This new distribution function, the g-Gaussian, was 
applied successfully to different astrophysical problems 
[7Hl6j. in addition of stellar radial and rotational veloc- 
ities. 

In the present work we study the behavior of the dis- 
tribution of stellar radial velocity in the solar neighbor- 
hood, using the Tsallis maximum entropy distribution. 
This work brings the analysis, for the first time, to our 
knowledge, of the behavior of the g-Gaussian distribution 
applied to a stellar parameter as a function of other ad- 
ditional relevant stellar parameters. In fact, we study the 
behavior of the g-Gaussian distribution of radial velocity 
as a function of stellar age, mass, rotation and distance 
of the stars to the Sun. This new investigation also of- 
fers the possibility of check the validity of the g-Central 
Limit Theorem, recently conjectured by Umarov, Tsallis 
and Steinberg [17] for physical sub-systems, such as those 
associated to stars and stellar systems. 

In section 2 we describe the sample and how we selected 
stars and physical stellar parameters for the present anal- 
yses. Section 3 brings the main results of the present anal- 
ysis. Finally, conclusions are presented in last section. 

Stellar Working Sample and Observational Data. 

— The present working sample is composed of 6781 single 
F and G spectral types dwarf stars, with precise radial 
velocity computed by Nordstrom et al. |18j . According 
these authors the sample is complete to the magnitude 
limit of 7.8 and should also be volume complete for the 
F- and G-type dwarfs stars to a distance of ^40 pc. The 
referred sample is defined in a volume around the Sun of 
^200 pc, as shown in fig. [T] 

For our analysis we taken the radial velocities, mass, 
age and projected rotational velocities from Holmberg et 
al. |19j . which gives an improved version of the work by 
Nordstrom et al. 18 , with a new calibration for all the 
relevant stellar parameters. All the selected objects are 
low-mass stars with effective temperature ranging from 
4500 K to 7800 K (T Q w 5778 presenting a mass 

range of O.65M to 2.0M Q with error estimation of about 
0.05Mq and metallicity range [Fe/H] between -0.75 and 
0.5. Nordstrom et al. [18] show that the distribution of 
their estimated metallicities obeys a Gaussian distribu- 
tion, with a mean of -0.14 and dispersion of 0.19 dex. For 
futher details on the observational procedure, calibration, 
and error analysis, the reader is referred to [T8lll9| . 

Results and Discussion. — In the following, we an- 
alyze the behavior of the distribution of the radial velocity 
in the context of the Tsallis statistics, taking into account, 
in particular, the dependence of such a behavior on differ- 
ent stellar parameters, namely age, mass, projected rota- 
tional velocity and distance of the stars to the Sun. The 
generalized g-Gaussian distribution function for radial ve- 
locity v r can be written as (e.g., [3j) 



Table 1: Best q and o" 9 -values determined using non-linear 
regression L-M method for all the stars of the present working 
sample. 



R 



0- 1 

1- 2 

2- 3 

3- 4 

4- 5 

5- 6 

6- 7 

7- 8 

8- 9 

9- 10 



386 
1701 
1939 
002 
452 
402 
325 
201 
22S 



32.78 
28.31 
31.24 
34.59 
32.22 
32.04 
33.37 
33.74 
34.74 
36.45 



1.37±0.37 
1.23±0.22 
1.26±0.21 
1.1±0.1S 
1.56±0.03 
0.75±0.24 
1.25±0.47 
1.46±0.0.39 
1.17±0.62 
1.65±0.75 



28.07 
25.87 
28.10 
33.39 
25.09 
44.98 
38.89 
32.76 
43.18 
35.95 



0.97 
0.99 
0.99 
0.99 
0.96 
0.97 
0.95 
0.97 
0.92 
0.01 



<!> g (v r ) = A q 



l-(l-g)- 



(4) 



where a q is the characteristic width and q is a free pa- 
rameter that denotes the entropic index. The parameter 
q is related to the size of the tail in the distribution (e.g., 
[16]). At the limit, when q — > 1, the above distribution 
reproduces the usual Gaussian distribution, given by 



(5) 



(v r ) = Aexp r - 



1 <3 = Sun 



All g-Gaussian distributions are two-parameters (g — a q ) 
nonlinear functions given by eq. (jj). The result of our 
computation for these two parameters, for well defined 
stellar age intervals is presented in tables Q] [2] and [3j 
which brings also the total number N t of stars in each 
age interval, UGauss of the Gaussian distribution from eq. 
([5]) and i? 2 -parameter obtained by a nonlinear regres- 
sion method based on the Levenberg-Marquardt algorithm 
20,21 . This method was used to compute the g-Gaussian 
with symmetric Tsallis distribution from eq. ([4]). The ob- 
tained values of the parameter g, strongly suggest that 
the distribution of radial velocity for the present sample 
of field dwarf stars is far from being in agreement with a 
standard Gaussian, since the values of this parameter are 
significantly different from 1, independent of the stellar 
age considered, as well as the considered spectral types. 
The error bars in fig. 2] correspond to a 0.05 confidence 
limit, as given in tables [2 [2] and [3] 

To illustrate the results listed in the referred tables, we 
present in figs. [2]and[3ja few exemples of the obtained fits 
for the distribution of stellar age and distance to the Sun, 
once we consider the usual and the g-Gaussian. As we 
can see, the g-Gaussian is the more suitable distribution 
fitting the obtained functions, irrespective of the region 
of the distributions, in contrast to the standard gaussian 
that fits very well only the central regions of the observed 
distributions. Such a behavior is the same, irrespective of 
the stellar parameter here considered. 

For instance, Carvalho et al. [3J have found that for 
stellar clusters older than 1 Gyr exists a positive gradient 
between the parameter q of radial velocities distributions 
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Fig. 2: Comparison between the semi-log plot of the observed (dots) and fitted Gaussian (dashed line) and (/-Gaussian (full 
line) distribution of radial velocity of F-type stars. Upper and lower panels represent, respectively, examples of the distributions 
segregated by age and distance to the Sun. 
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Fig. 3: Comparison between the semi-log plot of the observed (dots) and fitted Gaussian (dashed line) and (/-Gaussian (full line) 
distribution of radial velocity of G-type stars. Upper and lower panels represent, respectively, examples of the distributions 
segregated by age and distance to the Sun. 
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Table 4: Best q and er^-values determined using non-linear regression L-M method for the F- and G-Type dwarf stars for other 
stellar parameters. 



Parameters 


F-Typc 


G-Typc 


Massa(MQ ) 


9 




R 2 


q 




R 2 








0.90-1.10 


1.33±0.31 


31.31 


0.98 


1.43±0.03 


30.16 


0.99 


1.10-1.30 


1.21±0.08 


27.52 


1.0 


1.48±0.61 


31.32 


. 9 1 


1 .30-1.50 


1.27±0.13 


26.95 


1.0 


0.66±0.59 


39.71 


0.94 


1 50—2 00 


1.29±0.33 


26.15 


0.97 


1.23±0.1 


33.38 


0.88 




1.28±0.21 


1.2±0.33 


WS ini(km s" 1 ) 


9 




R 2 


<7 




R 2 


1-5 


1.2±0.21 


34.55 


0.99 


1.29±0.13 


34.95 


0.99 


5-10 


1.22±0.31 


31.44 


0.99 


1.45±0.45 


26.13 


0.96 


10-20 


1.13±0.20 


27.5 


0.99 


0.59±1.37 


34.42 


0.87 


20-40 


1.28±0.23 


24.74 


0.99 


1.62±0.30 


24.71 


0.87 


>40 


1.15±0.18 


26.36 


0.99 








average 


1.20±0.23 


1.24±0.6 


R*, (pc) 


9 


°q 


R 2 


q 




R 2 


0-30 


1.29±0.67 


24.44 


0.89 


2.0±0.32 


20.72 


0.99 


30-50 


1.29±0.30 


29.11 


0.96 


1.23±0.42 


34.80 


0.96 


50-70 


1.46±0.15 


24.56 


0.99 


1.26±0.33 


36.42 


0.97 


70-90 


1.34±0.17 


25.92 


0.99 


1.26±0.01 


33.50 


0.97 


90-110 


l.OlzbO.01 


33.91 


1.00 


1.42±0.67 


29.36 


. 1 


110-130 


1.40±0.28 


24.65 


0.98 


1.30±0.71 


30.59 


0.93 


>130 


1.33±0.27 


25.73 


0.98 


0.80±2.21 


36.70 


0.65 




1.30±0.26 


1.32±0.67 



Table 2: Best q and Uq-values determined using non-linear 
regression L-M method for the F-Type stars of the present 
working sample. 



Ago(Gyr) 



n- 



n-i 

1- 2 

2- 3 

3- 4 

4- 5 

5- 6 
>6 



200 
1573 
1571 
560 
217 
146 
156 



29.82 
27.96 
30.78 
33.75 
33.01 
43.21 
30.20 



1.43±0.45 
1.35±0.27 
1.26±0.21 
l.liO.OS 
1.57±0.54 
0.89±0.67 
1.29±0.33 



24.72 
24.19 
27.72 
32.50 
25.48 
44.82 
40.04 



0.92 
0.97 
0.98 
0.97 
0.88 
0.78 
0.97 



and the stellar age, while for young clusters (<1 Gyr) this 
correlation is small. From fig. |4] one observes that the 
values of q do not change whether we considered the whole 
of the stellar sample or the stars segregated by spectral 
types. In general, for the three cases shown in fig. [4J the 
values of the q are roughly constant over the stellar age 
and are consistent with the average value 1.28. Such a 
behavior for q is very similar to the results obtained by 
Carvalho et al. [3] for young stellar clusters. Interestingly, 
the great majority of the young clusters analyzed by those 
authors are located in the solar neighborhood near the 
galactic disc, namely at the same region where are located 
the stars of the present sample. 

In addition, we have analyzed the behavior of the q- 
parameter, obtained from de distribution of the radial ve- 
locity, as a function of stellar mass, projected rotational 
velocity and distance of the stars to the Sun. Table|4]sum- 
marizes the results from such an analysis, from where one 
can observe that the values of q are systematically differ- 
ent of 1, irrespective of the stellar parameter considered, in 
reinforcing that the Tsallis function fit the observed stellar 
residual velocity better than the usual Gaussian does. 



Table 3: Best q and cr 9 -values determined using non-linear 
regression L-M method for the G-Type stars of the present 
working sample. 



Agc(Gyr) 



R 



0- 1 

1- 2 

2- 3 

3- 4 

4- 5 

5- 6 

6- 7 

7- 8 

8- 9 

9- 10 



368 
342 
235 
256 
247 
240 
210 
146 



34.65 
30.38 
33.24 
34.31 
33.32 
38.91 
46.1 
31.34 
34.17 
36.97 



1.51±0.66 
0.93±0.13 
1.18±0.43 
1.37±0.03 
1.43±0.43 
0.84±0.71 
1.16±0.47 
1.49±0.19 
1.24±0.66 
1.78±0.83 



27.97 
31.0 
31.04 
29.46 
27.98 
41.23 
43.35 
30.17 
40.67 
34.12 



0.93 
0.87 
0.96 
0.94 
0.97 
0.92 
0.93 
0.96 
0.91 
0.89 



Conclusions. — We used g-Gaussian distributions to 
investigate the radial velocity of a sample of 6781 single 
field dwarf stars of spectral types F and G in the solar 
neighborhood. These low-mass stars have ages ranging 
between 1 and 10 Gyr and distance to the Sun from a few 
to tens of parsecs. In addition to the age and distance 
to the Sun, we studied the behavior of the distribution 
of the radial velocity as a function of two other relevant 
parameters, namely stellar mass and projected rotational 
velocity. The present analysis shows that radial velocity 
distribution for these stars is far from being a standard 
Gaussian model. The q-Gaussian distribution described 
in Tsallis nonextensive statistics is clearly the best distri- 
bution fitting the obtained functions. 

The relationship observed between the g-parameter and 
stellar age, stellar mass, rotational velocity and distance 
to the Sun indicates that the significant radial velocity 
deviation of the stellar sample from a standard Gaussian is 
a result of statistical dependence between the distributions 
associated with these variables. On the other hand, this 
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Fig. 4: Values of the fitted parameter q as a function of the 
stellar age for a) all stars, b) only F-Type stars and c) only 
G-Type stars of the present working sample. 

dependence reveals the effects of long-range interactions 
and the formation of high-energy tails consistent with the 
q-CLT, where the nonextensive character is observed. 
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